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Ab initio calculations predict a metal-insulator transition at zero temperature to occur in 
La4Ni308 at moderate pressures as a result of a pressure-induced spin-state transition. The spin- 
state transition that is seen at 105 K at ambient pressure from a low-temperature high-spin state to 
a high-temperature low-spin state has been observed to be shifted to lower temperatures as pressure 
is applied. From our calculations we find that a smaller unit cell volume favors the metallic low- 
spin state, which becomes more stable at 5 GPa. Similar physics should take place in the related 
compound La3Ni206, but on a different energy scale, which may account for why the transition has 
not been observed in this material. 



I. INTRODUCTION 



Spin-state transitions are observed in a vari- 
ety of systems as a result of the competition be- 
tween Hund's rule coupling {^h) and crystal field 
strength (Ac/). Typical examples are LaCoOa, 
where Co^+:d^ cations in an octahedral environment 
can occur in a non-magnetic low-spin (LS) state and 
various excited magnetic spin-states (both interme- 
diate [IS] and high spin [HS] have been predicted^) 
due to of the cation being the same order of 
magnitude as Ac/. Various Fe compounds also 
show spin-state transitions, and this is a common 
feature in organomet allies literature,^ where metal 
cations can be tuned to be in different spin states. 
It is common that the LS state is more stable at 
lower temperatures, but in some systems, the op- 
posite occurs.^ Important competition between Ac/ 
and Jh at high pressure has also been found in 
calculation of the Mott transition in MnO, where 
the insulator-metal transition, moment collapse, and 
volume collapse are found (experimentally and theo- 
retically) just above 100 GPa. The pressure at which 
the transition occurs is sensitive to competition be- 
tween these two energy scales, and not to the corre- 
lation strength to bandwidth ratio. ^ 

La4Ni308 is a low- valence nickelate that has re- 
cently drawn some attention^ because of the simi- 
larities in its crystallographic and electronic struc- 
ture with superconducting cuprates. It crystal- 
lizes in a layered structure^ formed by three Ni02 
neighboring planes, with these trilayers separated by 
fiuorite La/02/La layers, providing a highly two- 
dimensional electronic structure. This compound 
is obtained from La4Ni30io (Ref. 9) by eliminat- 
ing the apical oxygens neighboring Ni cations, thus 
the environment changes from octahedral NiOe in 
metallic La4Ni30io to square planar Ni04 in in- 
sulating La4Ni308. Two different Ni sites exist in 
the structure, depending on whether Ni is in the 



inner layer or the outer layers of the Ni02 trilay- 
ers. The compound undergoes a phase transition 
around 105 K that has been described via ^^^La nu- 
clear magnetic resonance measurements as a transi- 
tion to an unconventional low-temperature antifer- 
romagnetic (AF) phase. ^ Recent structural studies^ 
suggest the transition is due to a spin-state transi- 
tion on the Ni cations. 

Two possible spin states can occur in La4Ni308. 
The (on average) Nii-33+.(^8.67 cations sit in a square 
planar environment that leads to a large splitting be- 
tween the dx2_y2 and d^^ bands caused by the ab- 
sence of apical oxygens in this low- valence nickelate. 
Ac/ within the e^ doublet can be then comparable 
to if the former is larger, a LS state develops 
and if the latter is larger, a HS state would be more 
stable. This distinction is depicted in Fig. [TJ This 
simple model of the electronic structure of the com- 
pound suggests that the HS state would have a larger 
in-plane lattice parameter due to the slightly larger 
occupation of the dx'^_y2 orbital and also a smaller 
Ni-Ni inter-plane distance due to the de-occupation 
of the antibonding orbital, as has been recently 
measured.^ Moreover, the HS state leads to an in- 
plane AF coupling being more stable, consistent 
with the observations of its magnetic properties. 

As can be seen in Fig. [TJ the HS and LS 
states would lead to drastically different proper- 
ties. The HS ion leads to an insulating solution^^ 
caused by the formation of dz^ molecular orbital 
states. These are bonding-antibonding split around 
the Fermi level, due to their strong a-bond along 
the c-axis. However, the LS state would have a 2/3- 
filled dx2_y2 band around the Fermi level, leading 
to a metallic result. Experimentally,^ a reduction in 
resistivity above the phase transition at 105 K has 
been observed, which is consistent with at least an 
admixture of some Ni atoms in a LS state occurring 
at high temperatures. If all Ni atoms were in a LS 
state, the system would become metallic. 
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FIG. 1: Two possible spin states in the trilayer Ni com- 
pound La4Ni308 are represented depending on the value 
of the splitting between the d^^i and d^2_y2 orbitals and 
the Hund's rule coupling strength. On the left side, we 
see the high-spin state (larger Hund's rule coupling) and 
the low-spin state on the right (larger crystal field split- 
ting). Three Ni atoms are represented to account for 
the strong bonding along the c-axis and the formation of 
molecular orbitals with d^2 parentage. The Fermi level 
is represented showing the insulating HS vs. metallic LS 
solution. 



La4Ni308 and La3Ni2 06 (similar, but containing 
Ni02 bilayers instead of trilayers) are compounds 
close to a metal-insulator transition, and can be 
driven metallic by oxygen doping maintaining the 
same structureJ^ This transition has been recently 
analyzed for the La3Ni2 07_(5 series. For the case 
of La4Ni308, it has been shown that even though 
the resistivity is quite high and insulating-like, the 
NMR data has a metallic-like contribution coming 
from spin-spin scattering with quasiparticles at the 
Fermi surface (a constant term in the 1/TiT vs. T 
data^ fit for temperatures above the transition at 105 
K). In this paper, we explore the evolution of both 
spin states and their relative stabilizations, with re- 
spect to the different computational and physical pa- 
rameters (values of U, LDA+U scheme, volume vari- 
ations) suggesting that a spin-state transition can 
also be attained, with the corresponding reduction 
in resistivity, at moderate pressures. 



II. COMPUTATIONAL DETAILS 

Our electronic structure calculations were per- 
formed within density functional theoryi^ii^ using 
the all-electron, full potential code wien2k^^ based 
on the augmented plane wave plus local orbital 



(APW+lo) basis set.-^^ The generalized gradient 
approximation^^ (GGA) was used for the structure 
optimizations at each volume (that include both op- 
timizations of the c/a ratio and the atomic posi- 
tions) to estimate the equation of state of the sys- 
tem. Pressure values were obtained utilizing the 
Murnaghan^^ and also the Birch-Murnaghan^^ equa- 
tions of state. To deal with strong correlation ef- 
fects we apply the LDA+U scheme^SiSl including an 
on-site repulsion U and Hund's coupling J for the 
Ni 3d states. Results presented below compare two 
possible LDA+U schemes: the so-called "fully lo- 
calized limit" (FLL) and the "around the mean 
field" (AMF) scheme. A description of the results 
obtained at different values of U (in a reasonably 
broad range 4.5-8.5 eV for the Ni cations) is given 
in the main text below. The value chosen for the 
on-site Hund's rule strength is J = 0.68 eV. 



III. RESULTS 

A. La4Ni308 

The electronic structure of the HS state has been 
described extensively elsewhere ^iSiii^ Here we present 
the band structures of the LS state. This state, as 
has been calculated before,-^ has an in-plane ferro- 
magnetic (FM) coupling due to the less than half- 
filled dx2_y2 in-plane orbital, and the coupling out of 
the plane is small and AF. In Fig. [2] we present the 
band structure for one of the spin channels, with the 
"fat-bands" highlighting the e^ orbitals of the outer 
Ni atoms. The electronic structure in this phase was 
sketched in the right panel of Fig. [TJ where the d^s 
bands in both spin channels are shown to be fully oc- 
cupied. The left panel of Fig. [2] shows the minority- 
spin Ni bands highlighted. These are split into 
the bonding, non-bonding and antibonding bands. 
As described in detail in Ref. [lo| the inner Ni atoms 
do not contribute to the non-bonding molecular 
state (in Fig. [2] we have highlighted only the outer 
Ni for simplicity). We can observe that the partly 
filled dx2_y2 bands cross the Fermi level leading to a 
metallic state (these are almost degenerate for both 
the outer and inner Ni atoms), their bandwidth be- 
ing more than 3.5 eV. In previous works where the 
band structure of the HS state is presented, ^^^^^ the 
dx2_y2 bands are narrower, their bandwidth reduced 
by the AF in-plane coupling which is not stable in 
the case of a LS state due to the different filling of 
the dx2_y2 band. Dispersions along k^ are negligible 
in this highly two-dimensional compound. Hence, 
only the two-dimensional tetragonal Brillouin zone 
is analyzed. 

We have studied the energetics of the two possi- 
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FIG. 2: Band structures of the LS state are shown, calculated with LDA+U (around the mean field) for U= 6.8 eV. 
On the left (right) panel, the d^2 {d^2_y2) bands are highlighted. 2/3-filled d^2_^2 crosses the Fermi level leading to 
a metallic solution. 
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FIG. 3: Evolution of the total energy difference between 
the HS and the LS states with U for the two types of 
LDA+U schemes we have used: FLL (fully localized 
limit) and AMF (around the mean field). Positive (neg- 
ative) energy differences indicate the HS (LS) state is 
more stable. A crossover is found for the AMF scheme 
at around U= 6.5 eV. 



ble spin states. Figure [3] shows our analysis of the 
total energies we have computed. We present the re- 
sults of the total energy difference between the two 



spin states as a function of U (in a reasonable range 
of values from 4.5 to 8.5 eV) for the two LDA+U 
schemes used (FLL and AMF). FLL predicts the HS 
state to be more stable for all values of U. Its stabil- 
ity increases as U is increased. The AMF LDA+U 
scheme, however, tends to favor the stabilization of 
LS states. This is observed in Fig. [3l with the 
AMF method it is possible to obtain similar stabi- 
lization energies for the two spin states, whereas in 
the FLL scheme, an unrealistically small value of U 
would be necessary. As we can see from our calcula- 
tions, it is difficult to describe the energetics of the 
different spin states in La4Ni308 purely ab initio, 
we need to make additional assumptions. From ex- 
periments, a spin-state transition is observed in this 
compound at 105 K, to a low-temperature stabiliza- 
tion of the HS state. ^ Assuming that at temperatures 
above the transition the LS state can be thermally 
accessed, we could infer that the order of magnitude 
of the energy difference between both spin states is 
about 105 K (~ 10 meV/Ni). Accepting this, one 
can use the AMF scheme with a U 7 eV as a cor- 
rect description of the energetics of the spin state 
transition for this particular compound. This is the 
value and scheme utilized for the band structures 
presented in Fig. O 

It is of special interest to study which spin state 
becomes favored when pressure is applied to the sys- 
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tern. For doing this, we utilize the AMF scheme with 
U= 6.8 eV, a value that reproduces the energetics of 
the transition correctly, as we have just described 
above. It is difficult to estimate a priori which spin 
state will be favored at high pressure. The ener- 
gies coming into play will be the relative position of 
the two e^ states (Ac/) and the bonding- antibonding 
splitting in the bands. Applying pressure re- 
duces the Ni-Ni interplanar distance, this enhances 
the bonding-antibonding splitting (further promot- 
ing a HS state), but also lowers the relative energy 
of the band (a stronger metal-metal bond along 
the c-axis produces a larger d^^ occupation to screen 
the repulsion). Within the plane, reduction of the 
Ni-Ni distance destabilizes the dx2_y2 orbital favor- 
ing the LS state. All in all, several energies need to 
be taken into account and thus it is not easy to an- 
ticipate what pressure effects will do to this system 
before doing the calculations. 
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FIG. 4: Evolution of the total energy difference between 
the HS and the LS states as a function of volume. LS be- 
comes more stable at lower volumes. Calculations were 
carried out with the AMF LDA+U scheme for U= 6.8 
eV (see text). 



Results from the actual total energy calculations 
are summarized in Fig. HJ These show a reduction 
in volume favors the LS state. The calculations were 
carried out using the same atomic positions as ob- 
served experimentally^ and the same c/a value for 
all volumes considered. We observe a very rapid 
variation in the total energy difference, suggesting a 
transition from a HS to a LS state occurs at rela- 
tively small pressures. As we have discussed above, 
the LS state would lead to a metallic solution, with 
a 2/3-ffiled dx2_y2 band crossing the Fermi level. 
Thus a metal-insulator transition is predicted from 
our calculations (or at least the accessibility of LS 
states, that would eventually lead to metallic behav- 
ior) at a modest volume reduction of only about 1% 
with respect to the experimental one. 



To calculate the pressure at which we predict the 
transition to take place, we have carried out a series 
of GGA calculations for the structural optimization 
at each volume. GGA is known to give good esti- 
mates of the lattice constants, and this is generally 
not improved by the LDA+U method. Our GGA 
calculations predict a unit cell volume 1% larger 
than the experimental data (the theoretical lattice 
parameters are a= 3.95 A, c= 26.42 A, in good 
agreement with the experimental^ a= 3.96 A, c= 
26.04 A). We have fit our total energies calculated 
for each volume to both a Murnaghan^^ and a Birch- 
Murnaghan^^ equation, to obtain the pressure vs. 
volume relationship. The crossover between HS and 
LS, with the corresponding pressure-induced metal- 
insulator transition takes place at about 5 GPa, the 
bulk modulus of the compound being 180 GPa. 



B. Comparison with La3Ni206 

In principle, the physics of the spin-state transi- 
tion is the same in this compound and in La3Ni206. 
However, the latter does not show experimentally 
a spin-state transition, or this happens in a 
smoother way with no signatures of a transition in 
susceptibility, resistivity and heat capacity. Our pic- 
ture would predict the same metal-insulator tran- 
sition to occur on a different pressure/temperature 
scale, but the lack of experimental input about the 
energetics of the two possible spin states prevents us 
from being able to estimate the pressure required for 
the transition. For the same reason that it is difficult 
to predict what direction pressure will drive the rela- 
tive stabilization of a HS vs. LS state, it is difficult to 
compare the situation in La3Ni2 06 with La4Ni308. 
If we compare the band structures presented with 
the same value of U and the same LDA+U scheme 
(FLL) in Refs. |lo| and UH for the HS state, one can 
extract from them the band structure parameters 
(td^2-d^2 being the hopping parameter for the a- 
bond between Ni cations along the c-axis and Ac/ 
the intra-e^ crystal field splitting), that are sketched 
in Fig. [51 td^^-d^^ ^ 0-3 eV, A^/ - 2.6 eV for 
La3Ni206, and td^2-d^2 ^ 0-2 eV, A^/ - 2.0 eV 
for La4Ni308, which is consistent with their similar 
structural parameters (3.96 A for the Ni-Ni in-plane 
distance for both materials and a slightly larger out- 
of-plane Ni-Ni distance for La4Ni308: 3.25 vs. 3.19 
A). The larger crystal field splitting in La3Ni206 
would favor the LS state with respect to the case of 
La4Ni308, however the larger bonding-antibonding 
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FIG. 5: Comparison of the band structure parameters obtained for La4Ni308 and La3Ni206 in Refs. [13 andfl^ from 
LDA+U calculations. Observe that a smaller crystal-field splitting and a larger bonding-antibonding splitting occur 
for La3Ni206 compared to La4Ni308. Three (two) Ni atoms are represented for the case of La4Ni308 (La3Ni206) to 
account for the bonding-antibonding splitting of the d^2 bands. 
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FIG. 6: Evolution of the total energy difference between 
the HS and the LS states with U for the two types of 
LDA+U schemes we have used: FLL (fully localized 
limit) and AMF (around the mean field) for La3Ni2 06. 
Positive (negative) energy differences indicate the HS 
(LS) state is more stable. Comparing with the results 
for La4Ni308, the LS state is closer in energy, being al- 
ways favored in the AMF scheme. 



splitting in La3Ni2 06 goes in the opposite direc- 
tion favoring the HS state also. Thus while similar 
physical properties are highly expected in both com- 



pounds, the exact energetics of a possible transition 
in La3Ni2 06 are difficult to predict purely ab initio. 
One big difference is the existence of a non-bonding 
band in La4Ni308 just 0.3 eV above the top of 
the occupied da,2_^2 band (see right panel of Fig. 
[5|), whereas in La3Ni206 the antibonding d^^ band 
is significantly higher in energy, but the top of the 
valence band has dj;2_y2 character in that case (see 
left panel of Fig. [5|). 

To obtain a more qualitative understanding of a 
possible spin-state transition in La3Ni2 06, we have 
computed the total energy difference between the HS 
and LS state for the same two LDA+U schemes we 
used for La4Ni308 with the same range of U values 
and also in the experimental structure. In prin- 
ciple, the slightly smaller d-occupation (nominally 
d^"^ vs. d^"^^ in La4Ni308) suggests a very simi- 
lar but slightly smaller value of U should be rea- 
sonable for this compound. The results are sum- 
marized in Fig. [6l where by comparison with Fig. 
[3] for La4Ni308, it can be seen that the LS state 
should be more easily accessible for La3Ni206, sug- 
gesting that the larger crystal field splitting plays 
the biggest role in determining the relative stabi- 
lization of the different spin states. Using the AMF 
scheme, the LS state is always more favored, and 
the stabilization energy of the HS state in FLL is 
smaller than in the case of La4Ni308. However, no 
experimental indications to date suggest a spin-state 
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transition takes place in La3Ni2 06: neither specific 
heat^ nor resistivity^^ show any anomahes at a par- 
ticular temperature. Still, more experiments are re- 
quired to clarify the different properties of La3Ni206 
compared to La4Ni308. 

IV. SUMMARY 

Our ab initio calculations for the compound 
La4Ni308 describe the energetics of the spin- 
state transition predicted for the system. A low- 
temperature HS state has been observed experimen- 
tally and is predicted to be the ground state by our 
calculations at room pressure. These suggest the 
metallic LS state would be stable at moderate pres- 
sures (~ 5 GPa). Our electronic structure model 
for the system predicts a pressure-induced metal- 
insulator transition would occur together with the 
spin-state transition. Also, our description of the 



electronic structure can account for the different 
spin-state properties in La3Ni2 06, the changes in 
lattice parameters in La4Ni308 at low temperature, 
and gives indications of their proximity to a metallic 
state (for which some experimental evidence already 
exists for La4Ni308). 
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